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Abstract 

The next-to-leading order (NLO) QCD calculation for the isolated photon 
and isolated photon plus jet photoproduction at the ep collider DESY HERA is 
presented. The predictions for the isolated photon with no restrictions imposed 
on the jet are compared with the previous ones obtained in the small cone 
approximation, and the differences are found to be below 2%. The theoretical 
uncertainties in the cross section of the photoproduction of the photon plus jet 
are discussed. A short comparison with the new preliminary HI data and with 
the NLO predictions of Fontannaz et al. is also presented. 



1 Introduction 

In the previous papers [1, 2] we have presented the next-to-leading order QCD calcu- 
lation (NLO) for the photoproduction of a photon with large transverse momentum, 
so called prompt photon, at the ep collider HERA. We have also shown the leading 
order (LO) predictions for the prompt photon and a jet production [2] ^. 

For comparison with experimental data, the final photon was considered in [1, 2] 
as an isolated one, i.e. the hadronic energy within some cone around the photon was 
restricted. We introduced the isolation in an approximated way with the assumption 
that the momenta of partons inside the isolation cone were almost collinear with the 
photon momentum (small cone approximation). 

In this paper we study in NLO the photoproduction of both the isolated (prompt) 
photon and isolated (prompt) photon plus a jet without the small cone approxima- 
tion. Since the numerical differences between the present calculation and the previous 
approximated one are found to be very small for the photoproduction of the isolated 
photon, we focus our attention mainly on the NLO predictions for the isolated photon 
observed in the final state together with a jet, see fig. 1. 

The photoproduction is a process in which a real or almost real photon is collided 
with another particle in order to produce some final state (for a recent reviews see [3, 
4, 5]). In an electron-proton scattering the photoproduction corresponds to the small 
virtuality of the exchanged photon, Q^, say < 1 GeV^ [6, 7, 8, 9, 10]. In such 
a case this almost real photon colliding with the proton may lead for instance to a 
production of jets or particles with large transverse momenta. 

The photoproduction of the prompt photon in the ep or jp collision have been stud- 
ied theoretically (for a non-isolated final photon) since 1969 [11]- [24]. Then the NLO 
predictions for the isolated photon have been presented by various authors [25, 1, 26, 2]. 
The photoproduction of the isolated photon and a jet at HERA have been calculated 
in NLO by Gordon [27] and Fontannaz, Guillet and Heinrich [28]. There is also a 
NLO calculation for the isolated photon and isolated photon plus a jet production in 

rough estimation of NLO predictions for the photon plus jet production was given in [1]. 




Figure 1: The ep jet X photoproduction. 
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the deep inelastic scattering (Q^ ^ 1 GeV^) [29], which is in many theoretical aspects 
close to the calculations for the corresponding photoproduction events. 

The photoproduction of isolated photons without and with additional jets has been 
measured at the HERA collider by the ZEUS [6, 7, 8, 9] and HI [30, 10] Collaborations. 
In general, the data are in reasonably agreement with NLO predictions. Nevertheless, 
none of theoretical predictions as well as Monte Carlo simulations give sufficiently good 
description of the ZEUS data [7, 8] (see also [31]) for the isolated photon at negative 
photon rapidities (rij). Some discrepancies between NLO [28] and data are also seen 
for the preliminary HI data [10], mainly at r).y > 0.6 and < 0.5. 

In the paper [9] the ZEUS Collaboration has considered in Monte Carlo simulations 
the intrinsic transverse momentum of partons in the proton in photoproduction events 
containing the prompt photon plus jet, and the effective intrinsic transverse momentum 
was found to be < kr >— 1.69 ± 0.18l'lo!2o GeV. On the other hand, no need for the 
intrinsic transverse momentum in the proton was reported in [28]. 

In our analysis the intrinsic transverse momentum is not included, since we think 
that the measured < kT > is not due to the transverse momenta of partons inside the 
proton, but it describes effectively higher order emissions or multiple interactions. 

The paper is organized as follows. First, in sec 2 the processes contributing in NLO 
and the general formula for the cross section are shortly presented. The isolation cuts 
and the jet definition are introduced in sec. 3. In sec. 4 the slicing of the three body 
phase space is discussed, and some other calculation details and inputs are given in 
sec. 5. The predictions for the ep e'-fX process are given in sec. 6. The results for 
the ep — > 67 jet X process, discussion of theoretical uncertainties, and comparison 
with the new preliminary HI data, as well as with the predictions of Fontannaz et 
al. [28] are presented in sec. 7. A short summary is given in sec. 8. 

2 Cross section 

In our previous papers [1, 2] we have discussed the NLO calculation of the cross section 
of the isolated photon production in the ep collision. The set of contributing diagrams 
was different than in NLO calculations of other authors [25, 26, 27, 28], since we 
assumed that parton densities in the photon were of order 0{a), while in the cited 
papers their were assumed to be of order 0{a/as). 

The parton densities inside the photon are proportional to the electromagnetic cou- 
pling, a, and to the logarithm of the energy scale /i at which the densities are probed, 
Indj!^/ Aqqjj) ^ [32, 33, 34]. This logarithm is formally proportional to l/asin"^) but it 
originates from the pure electromagnetic process, 7 — > gg, with no strong interactions, 
and therefore the parton densities in the photon should be treated as quantities of 
order C(q;), and not 0{a/as). This different treatment of parton densities leads to 
different set of diagrams in NLO calculations involving resolved photons, and to some 
moderate differences in numerical predictions in comparison with predictions of other 
authors [2]. 

^We define the scale /x in sec. 5. 
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In the present calculation we take into account the same set of contributing partonic 
processes as it was taken in [2], namely: the Born diagram (7^' — > 75) together with 
QCD corrections, the box diagram {'yg — > ^g) [35], processes with the resolved initial 
or final photon, and processes with resolved both the initial and final photon. The 
general formula for the differential cross section for the ep — > e^X (or ep — > 67 jet X) 
photoproduction is: 

and 

a=l,q,q,9 b=q,q,g c=^,q,q,g'' ^ 

b=q,q,g 

where ah cd\ and 76 ^d\d2 are the partonic processes. 

The eq. 1 is the Weizsackcr- Williams approximation [36, 37] (see also e.g. [38, 39, 5]) 
with the flux of the real photons emitted from the electron given by G^/e(|/), where y 
is the fraction of the initial electron momentum carried by the photon. The functions 
fahi fb/p, s-nd D-y/c stand for parton distributions in the photon and proton, and a 
fragmentation function into the photon, respectively. For the direct initial (a = 7) or 
final (c = 7) photon the corresponding distributions arc replaced by the 5-functions: 
/a/7 = ^(2^7 ~ 1) OT Dry/c = 6{z — 1). The longitudinal-momentum fractions in the 
photon and proton parton densities, and in the fragmentation functions are denoted 
and z, respectively, while /i stands for the factorization/renormalization scale. 

The calculation of the partonic cross sections, da'^'^^'^'^'^^'^^\ involve isolation restric- 
tions and other kinematic cuts, as used in corresponding experimental analyses. 



3 Isolation cut and jet definition 

The NLO calculation includes partonic processes with two (2 2) or three (2 — > 3) 
particles in the final state: 

ab cdi{d2), (3) 

where a is the photon or a parton originating from the photon, 6 is a parton from the 
proton, c stand for the final photon or a parton which decay into the photon in the 
fragmentation process, and di are quarks and/or gluons. The outgoing partons di are 
not observed in the final state - they recombine into colorless hadronic jets. 

In experimental analyses the final photon is required to be isolated, i.e the trans- 
verse energy of hadrons inside a cone of radius R around the photon is assumed to be 
less than e times the photon transverse energy, 

El^'"^""'" < eE^, (4) 

hadrons 
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where e is a small parameter. The cone is defined in the rapidity and azimuthal-angle 
phase space, 



iVhadron " Vf^ + {(l^hadron - 0^)^ < R- (5) 

The some isolation restriction is taken into account in theoretical parton-level pre- 
dictions, where the sum in cq. (4) runs over c-parton remnant and over dj-partons, if 
they are inside the cone (5). In the present calculation we take R—1 and e = 0.1, as 
in the HI and ZEUS Collaborations' analyses. 

In the partonic 2 — > 2 processes with a direct final photon (0 = 7), 

ab — > 7^1, (6) 

the photon is isolated by definition. If the final photon comes from the fragmentation 
process (c 7^ 7), 

ab — > cdi, (7) 

it is isolated for the momentum fraction z > 1/(1 + e), and this cut is used in an 
integration in eq. (2). An inclusion of the isolation in the partonic cross sections for 
2 — > 3 processes, 

76 -fdid2, (8) 

is more complicated, since one needs to restrict the final momenta of two partons 
(eq. 4, 5), and to care about a cancellation and factorization of singularities, as de- 
scribed in the next section. 

For 2-^2 processes only one jet appears, and it originates from the rfi-parton. 
The jet's transverse energy, rapidity, and azimuthal angle are assumed equal to the 
transverse energy, rapidity, and azimuthal angle of the di-parton. 

Two partons, di and ^2, produced in the 2 — > 3 subprocess may lead to two separate 
jets or they may form one jet. A number of jets in the final state depends on the jet 
definition. In this paper the inclusive /cy-jet finding algorithm [40] is employed. 

Since in NLO calculation we deal with no more than two partons forming a jet/jets, 
the algorithm becomes very simple. If the distance between the partons, R12, defined 
as 

R12 = ^(?7<ii -r]'i^y + (0<^i -(f)d2)2 , (9) 

is larger than an arbitrary parameter Rj, then two separate jets arise with transverse 
energies, rapidities, and azimuthal angles of the dj-partons: 

E^t^^' = Et' , 7]^^'* =7]'^^ , (t)^"'^ = (1)'^^ , 1 = 1,2. (10) 
For i?i2 < Rj the rfj-partons are treated as ingredients of one jet with 

E't^' = E^' + (11) 

rf^i = (£;*r^<^i + E^^T]'^^) IE^\ (12) 

= (£;^i(/)<^i + £;^V'^')/^f *■ (13) 
Following experimental analyses, i?j = 1 will be used in numerical calculations. 
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4 Phase-space integration 



To obtain predictions, one need to perform an integration over four-momentum of at 
least one final particle, since unintegrated partonic cross sections contain 5*^^^ -functions. 
In case of 2 ^ 2 processes with no virtual gluon exchange this integration is straight- 
forward, and will not be discussed. 

There are three types of diagrams giving NLO corrections to the Born process, 

namely: 75 — > ^qg (real gluon emission), 751 — > ^qq, and the process "yq-^^q with 
a virtual gluon exchange. The cross sections of the above processes contains infrared 
singularities. For numerical integration over momenta of final particles, one needs to 
isolate these singularities [26, 27]. Like in earlier calculations for the ep —>■ e'y{jet)X 
reaction [26, 27, 28], also in the present one the phase space is sliced into a few parts. 

Let us assume that 9ji is an angle between the momentum of di parton and the 
momentum of an initial particle j, where j = e, p. Let us also define the distance, 
Ryi, between the parton and the final photon: 

R,^ = ^Jivd. + (0^. - (i>,r- (14) 

We define the parts of the phase space in the following way: 

• Part 1. In the first part the variable w (see Appendix A) is assumed in the range 
Wcut ^ w < 1, where Wcut is an arbitrary parameter close to 1: < 1 — Wcut <^ 1- 
This region of the phase space contains all types of NLO corrections: virtual gluon 
exchange, real gluon emission, and the process '-/g '-/qq. The virtual gluon exchange 
is a 2 — > 2 process, while the other ones are of 2 — > 3. However, for w close to 1, 
Wcut < w < 1, the two final partons in the 2 — > 3 processes are almost coUinear or/and 
one of the final partons is soft. For Wcut sufficiently close to 1 the kinematics of the 
2 — > 3 processes is almost the same as in the 2 — > 2 case: 

jb —>■ ■yd , d = di + d2, (15) 

where the "particle" d has four-momentum equal to the sum of four-momenta of the 
c/j-partons, and its mass is almost zero. In this case the final photon is isolated and 
the jet can be identified with the d- "particle" . 

The other parts of the phase space, parts 2-5 described below, are defined for 
w < Wcut, and contain only 2 — > 3 processes. 

• Part 2. In this region w < w^ut and 9cut > rnin{9ei^9e2)-, where 9cut is a small 
arbitrary parameter, 9cut ^ 1- Here, one of the final partons has the momentum 
almost coUinear to the momentum of the initial electron and it, for sufficiently small 
9cut-i does not enter the isolation cone around the final photon, and does not enter the 
cone defining the jet. The second parton has a large transverse momentum balancing 
the photon transverse momentum, so the final photon is isolated, and the jet consists 
(on the partonic level) of only this very parton. 

• Part 3. Here w < Wcut and 9cut > min{9pi,9p2). One of the final partons has 
momentum almost coUinear to the momentum of the proton, and, like in part 2, the 
final photon is isolated and the high-E'^ jet originates from the second parton alone. 



5 



• Part 4. Here w < Wcut and Rcut > min{R^/i, R.^^), where Rcut is a small parameter, 
Rcut -C 1. In this case one of the final partons, say di-parton, is almost coUinear to 
the final photon and the photon is isolated if < eEj,. This di-parton does not 
contribute to the jet. 

• Part 5. The last part is defined as the region in which there are no collinear 
configurations: w < Wcut, dcut < min{9ei,9e2), O^ut < rnin{9pi,9p2), and R^ut < 
'min{Rji, R^2)- 

In numerical calculations we apply in part 1 the formulae integrated over four- 
momenta of the final partons and virtual gluons [17, 41], namely the eqs. (24) and 
(37) from ref. [17]. In these formulae all the soft gluon singularities present in the 
virtual gluon and real gluon corrections are canceled, and all the collinear singularities 
are factorized into the parton densities. 

Parts 2-5 contain no contribution from the virtual gluon exchange, and no soft 
gluon singularities. Analytical expressions for cross sections corresponding to parts 2, 
3, and 4 are given in the Appendix B. All the collinear singularities appearing in the 
calculations are factorized into parton densities in the photon (part 2), proton (part 
3), and into fragmentation functions (part 4). 

The cross section in part 5 has no singularities at all, and one can perform an 
exact numerical integration over final four-momenta in any kinematic range, including 
isolation restrictions and other cuts. For this purpose we use the squared matrix 
element given in eqs. (36) (with e = 0) in ref. [17]. 

We assume that the cut-off angle 9c.ut is defined in the centre of mass of the initial 
photon and the initial parton. The predictions should not depend on a choice of 
unphysical cut-off parameters, 1 — Wcut) 9 cut and Rcut, if they are small enough. On 
the other hand, they can not be too small, since very low values lead to large numerical 
errors. The results presented in sees. 6 and 7 arc obtained with 1 — Wcut = ^cut = 
Rcut = 0.01. We have checked that the change of the predictions due to the variation 
of these parameters is neghgible, below 1%, if they are taken in the range 10~^ < 
1 - Wcut < 0.03, 10-^ < 9^t < 0.05 and 10"^ < R^at < 0.05. 



5 Other calculation details 

Our NLO calculations are performed in the MS scheme. The factorization/renor- 
malization scales in parton densities and fragmentation functions are assumed equal 
to the renormalization scale in the strong coupling constant and are denoted as /x. 
As a reference ^ equal to the transverse momentum (transverse energy) of the final 
photon is taken, /i — pj, — Ej,. For comparison /i — Ej^/2 and /i = 2£'j. will be also 
considered. 

The two-loop strong couphng constant, 

a ^ n _ 6(153 - 19iVj) ln[ln(/iVA^)] 

' (33 - 2iV^) ln(/x2/A2) ^ {?,?,- 2Nff ln(/i2/A2) J' ^ > 

is used with the QCD parameter A=0.386, 0.332 and 0.230 GeV for the number of 
active massless fiavours Nf=3, 4 and 5, respectively. In LO we use one-loop as with 
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A=0.123 GeV and Nj=A. The above A values are obtained by us from the world 
average of as at the scale Mz [42] : 

a^{Mz) - 0.1183 ± 0.0027. (17) 

To minimize theoretical and experimental uncertainties, o[s{Mz) was determined from 

precise data based on NNLO analyses only. The data given at scales different than 
Mz were extrapolated to the Mz scale using the four-loop coupling. Although we use 
one- and two-loop expressions, we have applied (17) as the best estimation of the true 
value of as{Mz)- 

In numerical calculations we take the Gliick-Reya-Vogt (GRV) LO and NLO par- 
ton densities in the proton [43] and photon [44], and the NLO GRV fragmentation 
functions [45]. For comparison we also use other parametrizations, namely Martin- 
Roberts-Stirling-Thorne (MRST2002) [46], CTEQ6 [47], Aurenche-Guillet-Fontannaz 
(AFG) [48] and (AFG02) [49], Gordon-Storrow (GS) [51], Cornet-Jankowski-Krawczyk- 
Lorca(CJKL) [52], Duke-Owens (DO) [15], and Bourhis-Fontannaz-Guillet (BFG) [53]. 

The initial electron and proton energies at the HERA collider are assumed equal 
to Ef. = 27.6 GeV and Ep = 920 GeV, respectively. In the Weizsacker- Williams 
approximation [36, 37] the photon spectrum in the electron is taken in the form [39]: 

with the maximal virtuality Q^ax = 1 GeV. The above formula describes the spectrum 
of equivalent real (transversally polarized) photons. We neglect longitudinally polar- 
ized photons and the interference between longitudinally and transversally polarized 
photons, since they give very small contribution [54]. As usual, we also do not take 
into account an emission of the final large-pr photon directly from the electron [55] . 

The results presented in next sections are obtained in NLO QCD with use of 
the GRV set of parametrizations with /i — Ej^ and Nf — 4 unless explicitly stated 
otherwise. 



6 Results for the ep e^X process 

The comparison with the ZEUS Collaboration data [8] for the photoproduction of 
isolated photons have been discussed in details in [2]. Recently a new preliminary 
data have been presented by the HI Collaboration [10]. For a comparison with these 
data, we apply kinematic cuts used in [10], namely the fraction of the electron energy 
transferred to the photon is restricted to the range 0.2 < y < 0.7, and the final 
photon rapidity and transverse energy are taken in the limits — 1 < ?7.y < 0.9 and 
5 < E} <10 GeV, respectively. 

In fig. 3 the differential cross sections da/dEj, and da/drj^, arc shown. Our exact 
predictions (solid lines) are compared with the predictions obtained in the small cone 
approximation (dashed lines) discussed in our previous paper [2]. The differences do 
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not exceed 2%, so the small cone approximation works very well, despite the fact that 
the isolation cone of radius i? = 1 is not a small one. 

The predictions are in good agreement with the HI preliminary data [10] (not 
shown) for both da/dEj- and da/df]^ distributions with exception of one experimental 
point at r]^ > 0.58, which lies slightly below the predictions. 

7 Results for the ep 67 jet X process 

There are two publications of the ZEUS Collaboration presenting results of measure- 
ments of the isolated photon plus jet photoproduction at the HERA Collider [6, 9]. 
In the first paper [6] the total cross section integrated over some kinematic range is 
given. The aim of the second one [9] was to study transverse momentum of partons 
in the proton, and no data for cross sections were presented (although some data for 
distributions of events, not corrected for the detector effects, were shown). 

In the new paper of the HI Collaboration [10] (see also [56]), the preliminary 
photoproduction data are presented for various differential cross sections of both ep — > 
e^X (considered above in sec. 6) and ep 67 jet X processes. 

As in sec. 6, wc impose kinematic limits used in the HI analysis [10]. The cross 
sections are integrated over 0.2 < y < 0.7, and —l<r]^< 0.9 and/or 5 < Ej^ < 
10 GeV with the jet rapidity and jet transverse energy in the range —1 < r]jet < 2.3 
and 4.5 GeV < respectively. If two jets are found with the above parameters, 
that with higher E'^f* is taken. 

7.1 Theoretical uncertainties 

As it is discussed in details in [28], the symmetric cuts for the photon and the jet 
transverse energy leads to unphysical results in next-to-leading or higher orders of 
calculations. This effect is shown in fig. 4a, where the LO and NLO predictions as 
a function of the photon transverse energy are shown in the Ej<-ia.nge wider than 
the range considered by the HI Collaboration. At Ej^ values close to the minimal 
jet transverse energy, -E^tmin = 4.5 GeV, the NLO differential cross section has a 
discontinuity: for (-EJ)- ~^ 4.5 GeV it has a strong maximum, and a minimum for 
(Ej,)^ 4.5 GeV. In the minimum the value of the cross section is even negative. 
This unphysical fluctuation is not present in the LO calculation. 

An integration of the differential cross section over the photon transverse energy 
higher than the minimal jet transverse energy, Ej^ > E^l^^^, leads to underestimated 
predictions in NLO. However numerically this effect is not very important. 

To avoid theoretical instabilities, one can consider a cross section averaged over 
some ii^y-bins, see flg. 4b. Note, that the NLO predictions are well deflned if one takes 
the cross section integrated/averaged over E]^ from E^^^^^^ — A to -E^f^j^ -|- A, if A is 
sufficiently large, say A > 0.3 GeV. So, the bins of a length 1 GeV presented in fig. 4b 
are large enough to avoid errors due to the symmetric cuts. 

The cross section for Ej, < E^^^^ is dominated by the NLO corrections, since the 
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LO contribution is suppressed by the E^r^ -cut (4.5 GeV < Ei^ ) and by the isolation 
requirement (sec. 3). We checked that the dependence on the renormahzation/factor- 
ization scale, ji, is not strong: variations of ji from Ej, to Ej^/2 or 2E^ lead to changes 
of the cross section less than 3% for Ej< < -Ejf^jj„ and up to 5% for E^ > E'^f^^^. 

The NLO predictions for various /i, and the LO predictions are also shown in fig. 5. 
In each presented here bin the dependence on the choice of is less than ±5% for both 
E't^' and rjjet distributions. The cross section is suppressed for E^* close to or higher 
than the maximal transverse energy of the photon, Ej, ,^^^^ — 10 GeV. At negative ?7jet 
the LO predictions are higher by 14%, and at positive r]jet they are lower than the 
NLO ones by 16-27%. 

Since for E'j. > Ej^^^^^ the difference between LO and NLO is not larger than 27% 
and the variation of NLO results due to the variation of the /i scale is small, up to 5%, 
the calculation is stable, and one can expect that the QCD corrections of higher orders 
are not sizable. On the other hand, for Ej^ < ^'y^mm the NLO corrections constitute 
almost 100% of the cross section, and the corrections of higher orders may change 
predictions. 

In fig. 6 the results obtained using different parton densities in the proton are 
shown. The predictions of CTEQ6 (NLO) [47] are 6% higher than the predictions 
of MRST2002 (NLO) [46]. The GRV (NLO) [43] densities give results higher than 
MRST2002 by 5-7% at negative rj^, and 3-5% at positive rj^. Differencies between 
CTEQ6 and GRV do not exceed 4%. 

The comparison between different parton densities in the photon is shown in fig. 7 
for the x°^^ distribution, where x°^^ is defined as 

x"''' = {E'j.^'e-'^'" + E].e">")/2yEe . (19) 

The GS (NLO) parametrization [51] give predictions lower than GRV (NLO) [44] by 
20-36% for < 0.9. This large difference is due to the charm threshold assumed 
in the GS parametrization at large scale, iJ? = 50 GeV^, much higher than thresholds 
in other considered herein parametrizations [2]. The results obtained with use of 
AFG (NLO) [48] and AFG02 (NLO) [49] are very similar, and only the latter is 
shown in fig. 7. It gives predictions up to 15% lower than GRV for x"'''^ < 0.9. At 
large-x°^* region, x"''''^ > 0.9, the cross section is dominated by the contribution of 
processes with direct initial photons, and the differences between predictions of various 
parametrizations are small. For the total cross section integrated over all x"''^ values 
the difference between GRV and AFG02 (GS) is 4% (16% ). 

We have also compared predictions of GRV (LO) [44] and the new CJKL (LO) [52] 
parton densities in the photon, as well as predictions of DO (LO) [15], GRV (NLO) [45] 
and BFG (NLO) [53] fragmentation functions. The total LO cross section (within the 
kinematic range considered by the HI Collaboration) for the GRV (LO) parametriza- 
tion is higher than for CJKL (LO) by 3%. The isolation requirement reduce the 
contribution of processes involving the resolved final photon [2] , so the dependence on 
the choice of fragmentation functions is weak, even if the fragmentation functions differ 
considerably from one another. The total cross sections obtained with DO and BFG 
(set I and set II) are lower than the predictions of GRV by 2% and 4%, respectively. 
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The GRV distributions for the proton, photon and fragmentation have been used as 
a reference, and each time only one parametrization has been changed for a comparison. 
The differences observed in the total cross section are not large (with an exception of 
the GS densities, which give predictions considerably lower than the other densities 
in the photon due to the specific treatment of the charm contribution). However, 
the differences can be larger if one changes simultaneously all the used distributions. 
For instance predictions of the MRST, AFG02 and BFG set are lower than the GRV 
predictions by about 11% . 

7.2 Comparing with HI preliminary data and with FGH pre- 
dictions 

The predictions shown in figs. 4-7 are in reasonable agreement with the preliminary 
data of the HI Collaboration [10, 56] (not shown), although some differences are 
present, especially for values below 6.7 GeV and for negative rjjet- 

For 5 < < 6.7 GcV and for rjjet < —0.3 the NLO result are higher than 
the measured cross section by two standard deviations. Smaller differences are also 
observed at positive r]^ values, where the predictions tend to overshot the data. In 
these kinematic ranges the predictions are in good agreement with the preliminary 
data, if one takes the number of active fiavours Nf = 3 instead of Nj = 4 ^, but 
the scale /i = Ej., for Ej, above 5 GeV, seems to be too large to neglect the charm 
contribution. 

Now, we make short comparison with the NLO predictions of Fontannaz, Guillet 
and Heinrich (FGH) [28] presented in [10, 56]. The FGH calculation takes into ac- 
count the QCD corrections to the resolved-photon processes, which are not included 
in the calculation presented herein (see sec. 2 and ref. [2]). They also use parton 
parametrizations different than the GRV ones used by us. 

In the considered kinematic range the total cross section of FGH is about 6% lower 
than our predictions. However differences are larger for the differential cross sections. 
For example, the FGH cross sections are lower by about 20% for da/dEj, and da/drj^'^^ 
at E^' < 6.7 GeV and t]-^"^* < —0.3, respectively, and they give better agreement with 
the data. On the other hand, our predictions for da/dr]^^^ {da/dx'^^) slightly better 
describes the data at rijet > 1-6 (0.25 < x"^^ < 0.5), where FGH predictions are 
about 45% (13%) higher than ours. In other kinematic ranges the differences between 
calculations are not larger than ±15% (usually below 10%), and both calculations lead 
to a similar description of the preliminary HI data, although the da/dE^r^* distribution 
of FGH is a bit more fiat and more consistent with data in shape. 

^The previous comparison with the ZEUS data [8] for the ep — > reaction led to opposite 
conclusions: the predictions for Nf = 3 were too low in each kinematic range and better agreement 
with data was observed for Nf = A and Nf = 5 [2] . 
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8 Summary 



We presented the results of the NLO calculation of the cross section for the photopro- 
duction of the prompt photon and prompt photon plus jet at the ep collider DESY 
HERA. 

A new exact calculation for the prompt photon production was performed. We 
found that the predictions agree within 2% with the results obtained in the small cone 
approximation used in our previous analysis [2]. 

The main aim of this paper was a detailed analysis of the prompt photon plus jet 
production in NLO. (In the previous analysis [2] we presented only the LO predictions.) 
The dependence on the choice of parton distribution functions was found to be of order 
10%. The uncertainties due to the variation of the renormalization/factorization scale 
are of order ±5%, and it may indicate that the corrections of higher orders are small. 
But for relatively low E'^, E'^ < E^^l^^^, the NLO corrections are very large while the 
LO contribution is suppressed, and we expect that in this region the corrections of 
higher orders could change the predictions. 

Our predictions are in reasonably agreement with the HI preliminary data [10, 56], 
nevertheless in some kinematic ranges there are discrepancies up to two standard 
deviations. 

Differences between our predictions and the predictions of Fontannaz, Guillet and 
Heinrich [28] presented in [10, 56] are usually below 10%, but in some kinematic limits 
they are larger. The largest differences, up to 45% at large rjjet, are seen for the da/dr]jet 
distribution. The experimental errors are too large to conclude which calculation gives 
better description of the HI preliminary data. 
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Appendix 



A The notation 




7(^7) 



9(P2) q(pi) 



Figure 2: An example of the 2^3 NLO process. The omission of the final gluon gives the 
lowest order Born diagram (^q ^q). 



Here we introduce the notation which is used in Appendix B and in sec. 4. There 
are two kinds of 2 ^ 3 partonic processes: 



and 



7(g) + q{p) j{p^) + q{pi) + g{p2) (see fig. 2), 



(20) 



(21) 



The four-momenta of particles contributing to the processes are given in brackets. Let 
us define variables v and w: 



1 + t/s , W = — M/(t + s), 



where 



{q + pf 



t = (q — p-yY and u = {p — p^Y . 



(22) 



(23) 



The V and w variables are in the range < f < 1 and < w < 1. Note, that for 
massless particles (pi + P2Y = sv{l — w), and in the limit {pi +^2)^ — ^ one obtains 
w ^ 1 {v can not be too low, since the final photon has large transverse momentum). 
For the 2^2 processes w = 1 hj definition. 
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B Cross sections of 2 — > 3 processes for collinear 
configurations 

In this Appendix we have collected analytical formulae for the QCD corrections to the 
Born process in the regions of the phase space labeled as parts 2, 3 and 4 (sec. 4). The 
formulae are derived with an assumption that parameter O^ut {Rcut) is small, < 9,.ut ^ 
1 (0 < Rcut ^ 1), and all terms of order 0{6^^f) are neglected for n > 1. 

All the collinear singularities are factorized into appropriate parton distributions with 
use of the dimensional regularization. 

The QCD corrections to the Born process have contributions from all parts of the 
phase space: 

= E E,^^ . (24) 

^ Pi \QCD corr i=l,2,3,4,5 ^ Pj \QCD corr 

As explained in sec. 4, in parts 1 and 5 we use formulae taken from ref. [17]. The 
corresponding cross sections for the 7p — > ■y{jet)X reaction in parts 2, 3 and 4 consists 
of the contributions of the ^yq — > ^yqg and ^yg — > ^yqq processes: 



E. 



7 \QCD corr. 



I d{s + t + u) V 
Jo ^JTl 



q=u,u.. 



fq/p{x)E^^- h fg/p{x)E. 



d^Pj d^p^ 



dx , (25) 



where i= 2, 3 or 4 and the summation runs over all quarks' and antiquarks' flavours. 
We include 2A^/ flavours in (^(72^3^''^^^, since there are Nf possible pairs qq, and in each 
pair the quark or antiquark can be collinear to the initial electron, initial proton or to 
the flnal photon. The expressions for the partonic cross sections in parts 2, 3 and 4 
are: 



E,^, = P,^,,(t/;,y£;e^,„OI^"(^)r , (26) 





1 


a 


d^Pj 




27r 


^^75-799- 


1 


a 


d^P-y 


vws^ 


27r 



E,^^ = -^-P,^,,{w,yEMM'^'^'»? , (27) 

fi^iq^iqg 1 rvc 1 - ?) 

E,-^-, = , ^ ^P,^,,{-^^,xE,e,^,)\M^'^-^'^{vw)\' , (28) 

d-^p^ [l — v)s^2'K ^ l — vw 

fi^ig^iqq 1 rvc ^ - v 



13 



[1 — V + vw)s^ 2Tr 1 — v + vw 



7 



1 ry — - 

P,_^,(l -v + vw, E^^R^t)\M''^-'^%- ■ )|2 , (31) 



— V + vwjs'^ zTT 1 — I) + 

where we have used the notation: 



^™-(^, E) = -\[z' + {l- zf] In ^ + 1 , (32) 



P^^,,-(z, £;) = 27Vc,eJP,^,,-(z, £;) , (33) 

P,_,,(^, £;) = In + 1 - ^} , (34) 

W.-)-?{^^ln^..), (35) 

|M^.-..-(^)|2 ^ aase f^.^'^'f , (36) 
^ ^(1 — 

Nc 

1^79-75 (^;)|2 ^ 2a2e^i^ , (38) 

|MT9^«^(^;)|2 = ^^|MT5-7g(i - ^;)|2 , (39) 
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